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Abstract. Let f be an arithmetical function. The matrix [f(i, j)] n xrt 
given by the value of f in greatest common divisor of (i,j), f((i, j)) as 
its i, j entry is called the greatest common divisor (GCD) matrix. We 
consider the generalization of this matrix where the elements are in the 
form f(i, 

1 Introduction 



The classical Smith determinant was introduced in 1875 by H. J. S. Smith p3] 
who also proved that 





(1,1) 


(1,2) •• 


• (1,n) 




det[(i,j)] nxn = 


(2,1) 


(2,2) •• 


• (2,n) 


= <p(1)<p(2) • • • <p(n 




(n,1) 


(n,2) •• 


• (n-,n) 





where (i, j) represents the greatest common divisor of i and j, and cp(n) denotes 

the Euler's totient function. 

The GCD matrix with respect to f is 



[f(U)]n> 



f((1,1)) f((1,2)) 
f((2,1)) f((2,2)) 

f((n,1)) f((u,2)) 



f((1,n)) 
f((2,n)) 

f((n,n)) 
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If we consider the GCD matrix [f (i, j)] n xn> where 

f(n)=£g(d), 

d|n 

H. J. Smith proved that 

det[f(i,j)] nxn = g(1)-g(2)---g(n). 

For g = cp 

f(i,j) = Y_ «P(d) = (M), 

d|(i,j) 

this formula reduces to ([1]). Many generalizations of Smith determinants have 
been presented in literature, see [T| [5l 171 [10| [13] . 

If we consider the GCD matrix [f(i, j)] n xn where f(n) = ^ d | n g(d) Polya and 
Szego [T2] proved that 

[f(i,j)]nxn = G-C T , (2) 

where G and C are lower triangular matrices given by 

g(j)> jl l 

0, otherwise 
and 

l ' 1 0, otherwise 
L. Carlitz [3] in 1960 gave a new form of ([2]) 

[f(i,j)]uxu = C n diag(g(1), g(2), . . . , g(u)) C£, (3) 

where C n = [Cij] nxu , 



D = [dij] nX n diagonal matrix 

di 



i, jU 

0, j /i ' 



gUJ, i- = ] 

0, i^j 



From ([3]) it follows that the value of the determinant is 

det[f(i,j)] nxn = g(1)g(2)---g(n). (4) 
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There are quite a few generalized forms of GCD matrices, which can be found 
in several references [2j O El [H [11] . 

In this paper we study matrices which have as variables the gratest common 
divisor and the indices: 



[f(i,j)]nxn 



f(l,(l,1)) f(1,(1,2)) ••• f(l,(l,n)) 
f(2,(2,1)) f(2(2,2)) ••• f(2,(2,u)) 

f(n,(n,1)) f(n,(n,2)) ••• f(n,(n,n)) 



2 Generalized GCD matrices 

Theorem 1 For a given arithmetical function g let 

f(n,m)=£g(d)- L 9(d). 

d|n d|(n,m) 

Then 

[f(i,j)lnxn = C n diag[g(1),g(2),...,g(n)]D^, 

where C n = [Cij] nxn , 

1, j I i 



^ "10, j /i ' 

d .._J 1» ) A 

dl) " \ o, j I i ■ 

Proof. After multiplication, the general element of 

A = [aijkxn = C diag[g(l ), g(2), . . . , g(n)]D T 

ag= ^ g(k) = ^g(d)- ^ g(d)=f(i,j) 
1^ | { d|n d|(n,m) 



IS 



Particular cases 

1. If g(n) = cp(n) then 



f(n,m) = J~ cp(d) - Y cp(d)=n 

d|n d|{n,m) 



[U, mj 



□ 
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We have the following decomposition: 



(t,j)lnxn 



1 

2- 
n 



(1,1) 
(2,1) 

(n,V 



1 

2 
n 



(1,2) 
(2,2) 

(n,2) 



1 

2- 
n 



d,n) 
(2,n) 



2. If g(n) = 1 then 

f (n, m) = t(tl) — t(tl, m) 

and 

[t(i) -T(i,j)] nxn = C n diag(l, 1,...,l) D n . 

3. Let g(n) = \i{n). From 

f 0, n = 1 

f(u,m) = ^~ — J~ (x(d) = < 0, n > 1,m > 1, (n, m) > 1 . 

d|n d|(n,m) [ — 1 , othewise 

we have 

[f (i, j)]nxn = C n diag(|0.(1 ), fi(2), . . . , (i(u)) D n . 

4. For g(n) = n, f(n, m) = cr(n) — a((n, m)) and 

[f(i-,j)]nxn = C n diag(l, 2,...,n) D n . 

Remarks 

1. Due to the fact that the first line of the matrix [f(i, j)] n xn contains only 
0-s, the determinant of the matrix will always be 0. 



2. We can determine the value of the matrix associated with f , if the function 
f is of the form 

f (n, m) = H(n) — H((n, m)) . 
By using the Mobius inversion formula, we get 

g(u) = 2>(d)hQ), 
d|n 

consequently by using Theorem [TJ the matrix can be decomposed according 
to the function H(n): 

[f(i, j)]nxn = C n diag[(u- * h)(l ), {[i * h)(2), h)(n)] D n . 



164 



A. Bege 



Theorem 2 For a given arithmetical function g let 

n 

m, j) = J_ g(k) - Y. 9(d) - L 9(d) + ^ g(d). 

k=1 d|i d|j d|(i,j) 

T/ien 

nxn — 

w/iere D n = [d i j] nxn , 

dl > " \ o, j 1 1 ■ 

Proof. After multiplication, the general element of the matrix 
A = [a ij ] nX n = D n diag[g(1 ), g(2), . . . , g(n)] D n 



is 



£ g(k)=£g(k)- ^ g(k) 

j n k=1 k|n or k|m 

k /m 



n 



^g(k)-^g(k)-^g(k)+ £ g(k) = f(i,j). 
k=1 k|n k|m k|(n,m) 



Particular cases 

1. If g(n) = <p(n) then 

n 

f(U) =^(p(k)-i-j + (i,j), 
k=1 

[f(i,j)] nxn = D n diag[(p(1),(p(2),...,cp(n)] D n . 

2. If g(n) = 1 then 

f(i,j) =n-T(i) +T(i,j) 

and 



□ 



[f(t,j)]nxn = D n diag(l, !,...,!) D n . 
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3. g(n) = n. Then 

_,. n(n + 1) 

= — 2 a ^ ~ °"^ n ' m ^ 

and 

[f(i-,j)]nxn = D n diag(l, 2, ...,n) D^. 
Another generalization is the following: 

Theorem 3 For a given arithmetical function g let 

n 

mo) = y_ gw - Y. 9 (d) - L 9 (d) + L 9 (d) - 

k=1 d|i d|j d|(i,j) 

We define the following A = [aij] nX n matrix 

f f (U)> i->j>1 
' I g(D + m,j), i = 1 or j = 1 ■ 

T/ien 

A = D; diag[g(l),g(2),...,g(n)] D; T , 

w/iere = [d[j] n><n , 

1, 1 = 3 = 1 



dij, ij / 1 



Proof. 

We calculate the general element of the matrix 

B = [cHjkxn = d; diag[g(l), g(2), . . . , g(n)] D; T . 
If i > 1 or j > 1 we have 

n 

^ j/ n k=1 k|n or k|m 

k }m 



If i = j = 1 



^ g (k)-^g(k)-^g(k) + £ g(k) = ay. 

k=1 k|n k|m k|(n,m) 

bn = g(1) = an. 
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□ 



Particular cases 

1. If g(n) = cp(n) then 



n 



cp(k) -i- j + 



i,j > 1 



n 



^<p(k)-i-i + (i,j) + 1, i 



1 or ) = 1 



k=1 



2. If g(n) = 1 then 




n-x(i) -t(j) + x(i,j), 
n-T(i) + x(i,j) + 1, i 



i,j > 1 
1 or j = 1 



The following problems remain open: 

Problem 1 Let F(n, m) 6e an arithmetical function with two vriables. Deter- 
mine the structure and the determinant of modified GCD matrices 
A = [a(i,j)] nxn , where 



Problem 2 Determine the structure and the determinant of modified GCD 
matrices A = [a(i, j)] n xn> where 
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